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1. INTROOUCTIUN 


1.1 BACICSRQUNO 

The objective of this paper Is to present the more common vegetative indices 
used with Landsat remotely sensed data and describe their origins, develop- 
ment, logic, and relationships to ground-based measurements of vegetation. 

An effort has been made to preserve the order In which the various vegetative 
Indices appeared In the literature In order to historically trace their 
underlying concepts and develo usrt. A brief discussion of remote sensing 
preprocessing techniques as Input to the vegetation Indices is also included. 

A vegetative Index Is simply a formula that transforms the four-dimensional 
Landsat data Into a single real number. This number, the vegetative index 
number (VIN), may be the reflective count of a single Landsat band, or more 
commonly a combination of reflective counts of 2 or more bands. In general, 
the greener and denser the vegetation Is in an area, the higher the VIN. The 
Idea has been to formulate VINs that can be used to predict general crop 
health, crop growth stages, and crop yield. 


Since VINs are computed using reflection measurements of one or more Landsat 
bands, some discussion of the Landsat bands follows. Landsat Is a multi - 
spectral remote sensing satellite which receives and records reflectance in 
four spectral bands or channels. The wavelengths and spectrum of these four 
bands are presented below: 


CHI » Band 4 
CH2 - Band 5 
CH3 « Band 6 
CH4 » Band 7 


Wavelength, gm 

0.5 to 0.6 
0.6 to 0.7 
0.7 to 0.8 
0.8 to 1.1 


Spectrum 

Visible green 
Visible red 
Near-Infrared 
Near-Infrared 
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Bands 4 and 5 are called the visible bands, and these tMO bands generally nave 
lOM reflectance values where there Is live green vegetation, p<!'rt1cu1ar1y 
band 5. Bands 6 and 7 are the Infrared bands, and have high reflectance values 
where there Is live green vegetation. Bauer (ref. 5, page 5) sumoarlzes the 
reasons for this very concisely: “Ihe 1(M reflectance and transmittance of 
visible radiation Is attributed to the high absorption by leaf pigments, 
primarily the chlorophylls. However, these pigments are highly transparent to 
Infrared "adlatlon, and the Internal cellular structure of the leaf appears to 
determine the high reflectance at the<^e wavelengths." As a matter of fact, It 
Is known that vegetation reflects more In the near>1nfrared bands than do most 
other natural objects (ref. 7, page 2-2), making the near-infrared bands 
apparent Imniadlate Indicators of areas of vegetation, percent ground cover, and 
biomass. Landsat satellites up to this time have not recorded reflectance 
below .5 un or above 1.1 pm, a fact which limits research in some important 
areas of stress (see section 3.2). 

Some further properties of both individual leaf reflectance and crop canopy 
reflectance which are useful for the Landsat user to know are below, and are 
taken directly from Bauer (ref. 5, page 5-6). For individual leafs, other 
important factors affecting reflectance are maturation, senescence, water con- 
tent, nutrient stress, disease, and insect infestation. In general, as leaves 
mature, reflectance decreases in the visible spectrum and increases in the 
infrared spectrum. On the other hand, senescence produces the exact opposite 
reflectance responses. Both visible and infrared spectral reflectance 
Increases as leaf water content decreases; however, changes in reflectance are 
not substantial until the leaves reach about 75 percent turgidity. Thus, the 
change In reflectance Is not a sensitive indicator for Initial stages of 
drought. Nutrient stress affects reflectance in both the visible and infrared 
wavelengths, but Increases or decreases in reflectance are dependent upon the 
type of nutrient stress (see section 3.2). Disease and Insect Infestation are 
also known to affect reflectance; however, the wavebands which are necessary 
for specific detection of these problems on a large scale are restrictive to 
the extent to keep this area mostly descriptive (see ref. 7, page 4-4; ref. 9, 
page 1160). 
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When changing from Individual leaf cpactral maasureflients to larger area 
remotely sensed crop canopies, new variables arise which must be considered; 

1. Variations In mount of leaf area, biomass, and $^nd cover (hje to 
differences In planting date, soil type, soil moisture, plant population, 
and/or disease conditions. 

2. Variations In maturity (hie to differences In variety, planting date, soil 
type, and soil moisture. 

3. Differences In cultural practices, such as tillage or harvesting. 

4. Geometric configuration of the crop due to differences In row width, row 
direction, or lodging of plants. 

5. Environmental variables, such as atmospheric conditions, wind, angle of 
reflection In relation to solar Incidence angle, and soil moisture 
conditions. 

Part II of this paper details the mathematical derivations of a number of VINs 
used today and Involved in this is some basic logic or VIN theory. In all 
derivations, a few basic facts provide nuch of the background reasoning which 
will be briefly summarized from Tappan (ref. 19, pages 18-20). The wavelength 
of .68 is a very good single wavelength to discriminate between living vegeta- 
tion and dead or dormant vegetation when using wavelengths between 0.4 and 
1.1 ym such as Landsat does. Generally, the dead or dormant vegetation has 
higher reflectance than the living vegetation In the visible portion of the 
spectrum and lower reflectance In the near-infrared portion. The .68 wave- 
length also appears to discriminate well between soil and living vegetation, 
although in this case the optimum wavelength has been found to vary from 
environment to environment. Living vegetation has minimum reflectance values 
in the .35 to .5 and .67 to .69 um ranges, and maximum reflectance values in 
the .8 to 1.1 range. Referring back to page 1 at the chart on Landsat bands. 
It Is seen that Channel 2 (or Band 5) contains the Important .67 - .69 ym 
range, and that Channel 4 (or Band 7) contains the .8 - 1.1 ym range. 
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A« ntntloncd tirtltr. nest VlNs uietl today Involve conOlnatlons of Undsat 
bands ratiMr than single Undsat bands. It Is often hoped that a particular 
combination of bands Mill yield more Information than can be discerned from 
Individual band values. Multi band VI Ns tend to be more stable, and thus pro« 
vide better capability for season to season comparisons of vegetation amounts 
and conditions (ref. 14, page 1550). Tappan (ref. 19. page 27) gives another 
important advantage of using multiband VINs rather than single band VIds: 

In many cases It reduces the amount of "noise" In the Undsat system. Factors 
such as changes In atmospheric conditions betMcen Images taken at different 
times of the year, attenuation effects due to clouds and haze, topography, 
shadOM, soil and dead vegetation and others may affect all four Undsat bands 
In either similar or dissimilar Mays, but In many cases a band ratio or 
normalized band coart>1nat1on H11 result In the partial reduction In the noise, 
soflwtimes almost complete reduction. 

In summary, VINs are used to clarify Information content of Landsat bands for 
a defined purpose. Purposes of Interest Include (1) estimates of biomass, 
leaf area Index, density of ground cover and plant height; (2) Identification 
of stressed areas and stress factors; and (3) general and specific crop 
identification . 

1.2 PREPROCESSING OF LANDSAT DATA 

It was mentioned In the Introduction that Landsat data is altered by sun angle 
differences, clouds and haze, shadow, and other factors. Some researchers 
have given attention to the Idea of preprocessing Landsat data to reduce the 
effects of these factors before making use of the data. As previously 
mentioned, a band ratio or normalized band con4>1 nation will often result in 
the partial or almost complete reduction In sun angle or noise factor effects 
without requiring preprocessing. Nevertheless, It Is obviously desirable to 
have data which contains as little noise as possible before being used In any 
sense. Some of the preprocessing techniques are given below. 
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Tht cotint fun tnglt orrtctloi tlgorlthm (rtf. 12, ptgt 710) Is a nathaMtlcs 
transformation ««(t1ch corrtcts Landsat data to a rtfartnct solar tltvatlon 
angit. Tht corrtc.lon Is appiltd as folloMs: 

Ltt ■ Landsat signal In band f, 9 ■ solar ztnith angit, 9^ ■ rtftrtnct 
solar ztnith angit, and Xj - tht corrtcttd Landsat band 1 signal. 

Then 


cos 9 


XI • 


1 coTT ''1 


All resulting : will appear to have been acquired at the reference solar 
ztnith angle. 


The XSi'AR haze correction procedure, developed by Lamoeck (ref. 12. page 71d; 
ref. 8, pages 40>43), Is a haze correction procedure which has been undergoing 
revision since 1977. Tht first version, now known as the global XSTAK proce- 
dure (ref. 12), Is an algorithm which can be easily applied to Landsat data. 
However, It has generally been replaced In favor of the new spatially-varying 
XSTAR procedure (ref. 8). The spatially-varying XSTAR procedure is a multi- 
step software procedure, and program products from it, along with information 
about the software, are available from ERIM (Ann Arbor, Michigan). 

J. Potter (ref. 15) developed the Atmosphere Correction (ATCJR) program in 
1977 to correct Landsat data for naze, sun angle, and background reflectance. 
The program logic Is based. In parts, on principles of radiation transfer 
theory. The program Is available for use at the LARS coti^iutatlon center (West 
Lafayette, Indiana) and at the NASA JSC Bldg. JI7 Computation Center (Houston, 
Texas). Information concerning the use of ATCOR may be obtained by contacting 
Lockheed Engineering and Managetwnt Services Co. (Houston, Texas). 


In suMMiry, It sNwId be noted thit none of the (tevelopment of the VINs pre- 
sented In this paper considered preprocessing procedures, jihI as a matter of 
fact, some VINs had noise leductlon Included In their development. Of course, 
afiy band ratio VIN Implicitly corrects for sun angle due to cancellation in 
the division (ref. 7, page 3-2). 
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2. OONCEPlUAL AND MATHEMATICAL OERIVATIONS OF illii 


2.1 PEARSON'S AMO MILLER’S MULTIBAMO ViHs 

Although single band values from Landsat can be considered VINs. many 
researchers refer to single band VINs as simply single band values and only 
regard nultiband VINs as true VINs. Ihe rest of this paper ^IM gi^nerally 
follOM this outlook. Ihe pioneers of VINs In this context are Pearson and 
Miller, In 1971, and since then many VINs have been formulated. 

Pearson and Miller conducted extensive studies In 197 1*72 on vegetation canopy 
reflectance using grass plots (ref. 14; ref. 19, page 21 for summary of their 
findings) and among other conclusions chose .6H um and .76 u«n as the two opti- 
mal wavelengths for separating green vegetatlwi from soil and dead or dormant 
vegetation. Figure 1 Is a copy from page 1362 of ref. 14 of a graph wnicn 
shows the reflectance curves of soil, dead or dormant vegetation, and live 
green vegetation superimposed for what was considered a typical snort grass 
prairie plot. They also found that an Inverse linear correlation existed 
between reflectance and green vegetation at .68 ixn, and a direct linear 
correlation existed between reflectance and green vegetation at .78 i^iu Uoing 
further, they found that using combinations of the reflectance of these two 
wavelengths yielded even higher correlations between reflectance and total 
biomass (live and dead biomass) _ sometimes r > 0.90. First, they correlated 
total biomass with .78 minus .68 un reflectance values and obtained an 
r - .88; an r ■ .91 was obtained when they correlated total biomass with the 
ratio of the reflectance at .78 um to that at .68 ixn. The underlying princi- 
ple In both cases Is that as the quantity of green vegetation Increases in a 
given area, the red wavelength (.6 to .74 t^) reflectance decreases whereas 
the near-infrared wavelength (.74 to 1.35 ivn) reflectance Increases. This 
same principle Is used In many of the subsequent VINs which quickly followed, 
many of which are still used today. 
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Figure 1.- Typical wavelength versus reflectance curves 
(Source; ref. 14, p. 1362). 

2.2 DEVELOPMENT OF THE TV I 

As part of the Great Plains Corrider Rangeland Project conducted at Texas A4M 
University, Rouse et al. (ref. 17) In 1973 developed a VIN which they called 
the Transformed Vegetation Index (TVI). It uses the basic theory developed 
only a year or two earlier by Pearson and Miller: that a contrast between the 
.78 and .68 pm wavelengths provided a better correlation with biomass than did 
single wavelengths. 

The work carried on In this project was some of the first to use Landsat data. 
Since Landsat does not have Individual spectral wavelengths (such as .68 pm 
and .78 pm) but rather wavebands. Rouse and others use bands 5 and 7 (or chan- 
nels 2 and 4) to correspond with the wavelengths .68 and .78 pm respectively. 
The contrast used was the spectral value of channel 4 minus the spectral value 
of channel 2. Sun angle and haze were seer nt be a problem, and a normaliza- 
tion procedure was used to eliminate the effects as much as possible: dividing 
the contrast of the two channels by the sum of the two channels. The addition 
of 0.5 to this quotient was done to avoid working with negative numbers. Since 
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thty thought that the variance of this result night oe prt^ortloned to the mean 
values, they took the square root of It. Ihe result was named the TVI. The 
above verbal Ascription of the computation of the TVI Is simplified as 
follows: 


TVI - 




dHa - CHZ ^ 
CH4 V m * 


0.9 


Further analysis of their data Indicated that the TVI was not sensitive to 
sparsely vegetated areas. However, results of the study also concluAd that 
the TVI was adequate for monitoring the vernal progression and retrogradatlon 
of vegetation within the Great Plains Corr<Ar, (ref. 17, page 313), ar.:i also 
had potential for measuring Increments In green blc^ss, useful In regional 
agricultural applications. 

The TVI Is often now referred to as the TVI 7. The reason for this Is tnat 
another TVI has been developed, the TV 1 6, which Is the original TVI using 
Landsat hand 6 (channel 3) In place of band 7 (channel 4). Since the 
1tr^)ortant .78 kh wavelength Is seen to be bordering bands b and 7, the 
.78 un • .68 un contrast can be attained using either band, althougn so*ne 
prefer band 7 over band 5 and some vice-versa. In summary, Tvl7 and TVIo are 
computed as 'ollows: 


TVI 7 
7V16 


rCTW"‘-CH2 "■ 
“K rw-TTH? ^ 

I CH3 - Cfl2 ■ 
■ 1 CH3 * CHZ ' 


U.5 


0.5 


2.3 OEVELOPMENT OF il€ AVI 


About the saiT* time that Rouse and others were conducting their studies wnicn 
led to the creation of the TVI, other research was being conducted by 
P. Ashburn (ref. 3) which led to the creation of the Ashburn Vegetation Index 
(AVI) In 1974. Ashburn's Intention was to provlA some measure In green grow- 
ing vegetation, and the hope was that this measure could be helpful in crop 
1 Antificatlon studies using ViMs; section 11 of this paper reports on the use 
of the AVI toward this goal. 
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Like the TVI, the AVI makes use of the basic theory developed i years earlier 
by Pearson and Miller: that the .78 urn - .68 pm wavelength contrast provides 
a better greenness measure than do single wavelengths. Using Landsat data, 
Ashburn set the AVI « 2(CH4) - CH2, where the doubling of channel 4 normalizes 
the two channels (since CH4 digitizes from 0-63 for black to white and all 
other channels digitize from 0-127 for black to white). All negative AVI 
values are set to zero. 

Two things about the AVI make Its use very desirable In large-scale computer 
applications. The first is the simplicity of the formula. A simple calcula- 
tion is very desirable when processing many Landsat segments _ each of which 
have 22,932 pixels with 4 channels of Information. The second is the dichoto- 
mous implication which results from AVI computation: a positive AVI signifies 

at least some growing vegetation in a scene, a zero AVI signifies no growing 
vegetation. This makes the AVI a favorite for masking applications. 

2.4 OEVELOPMENT OF THE GVl, S8I, KVI, AND GIN 

In 1976, R. J. Kauth and G. S. Thomas published a paper (ref. 13) the results 
of which have had great impact on agricultural remote sensing research. In 
this paper, they outline a transformation which can be applied to Landsat data 
which preserves four independent dimensions of a Landsat scene. The remainder 
of this section is devoted to describing the development of this transform and 
the VINs which result. 

Their work began by inspecting scatter plots of 1973 channel 2 versus chan- 
nel 3 digital values obtained from Landsat 1 over certain Illinois counties 
(ref. 12, page 705), The scatter plot's values were not the digital values 
for single Landsat pixels over a segment, but instead were approximate digital 
values of numerous clusters of Landsat pixels for several segments. For each 
of these clusters, the exact soil type, vegetation cover, etc. for the actual 
land areas in each cluster were known. 
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A shape resenbllng Mhat Kauth and Thomas called a tasselled cap was seen to be 
ttw overall pattern In the scatter plots. Figure Z displays an actual scatter 
plot In Illinois and an outlined tasselled cap. It was seen to be a typical 
scatter plot and used for further study. The clusters along the base of the, 
cap In figure 2 were known by Kauth and Thomas to be varying shades of bare 
soil. Thus, the base of the cap could be used as a soil brightness line. 

This soil brightness line Is, of course, seen only In two-dimensional space 
since figure 2 Is two-dimensional (channel 2 versus channel 3). The next 
question Is what happens when three dimensions are considered Instead of just 
two? Would the soil brightness line change Into a soil brightness plane with 
an additional Landsatl channel considered, and what additional changes with 
the fourth channel also considered? When Kauth and Thomas added channel 1 to 
their scatter plot study of the figure 2 clusters, they found a shape which 
resembled a tasselled cap In three-dimensional space (this tasselled cap 
appeared to be one which had not been opened to wear — it was still basically 
planar but not completely). Figure 3 is a resemblance of the three- 
dimensional (channel 1 versus channel 2 versus channel 3) tasselled cap. 
Because of the small "thickness" channel 1 added to the tasselled cap 
representation, a thin cigar-shaped plane of soils resulted in the three- 
dimensional plots; this resulted In Kauth and Thomas continue to think in 
terms of a basic line of soils, even in the three-dimensional space (see 
fig, 3). The last channel, channel 4, was found to have an almost identical 
graphical effect to that of Channel 1. 

With the various scatter plots described above, Kauth and Thomas w_>re now able 
to describe a line of soils In four-channel space,, This is done by connecting 
the point (0, 0, 0, 0) to the point Rj (In Flgurf* 3), where Rj is obtained, of 
course, from looking at the numerous scatter plots that had been created, 

(even though Rj in fig. 3 Is only shown in its three-dimensional space). The 
description of this line of soils is simply a vector, or soil brightness 
vector. As a unit vector, the soil brightness vector is described as follows: 

.433(CH1) + .633(CH2) + .586(CH3) + .264(CH4) 
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Figure 2.- Cluster patterns from Fayette County, Illinois, 
June 11, 1973 (Source: ref. 13, p. 4B-47). 


Channel 3 



Figure 3.- Three-dimensional Tassel led Cap representation 
(Source: ref. 12, p. 707). 
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For aoy pixel, a soil brightness measure could be (H)ta1ned by using the above 
soil brightness unit vector with the CHI - CH4 digital values of that pixel. 
This soil brightness measure Is iww called the Soil Brightness Index (SBI), 
and as seen from above, 

SBI - .433(CH1) + .633(CH2) + .586(CH3) + .264{CH4) 

At this point, Kauth and Thomas reexamined the tassel led cap scatter plots and 
noted that the clusters which contained pixels having high greenness or living 
vegetation appeared toward the top of the tassel led cap. The farther the 
cluster was from the line of soils of figure 3, the more vegetation or 
greenness In that cluster. Thus, they created a new vector pointing in a 
perpendicular direction from the line of soils to R 2 ; this vector was actually 
generated by using the Gram-Schmidt orthogonal izat ion procedure and was also 
made into a unit vector. This vector, which is obviously considered to be a 
measure of greenness, is described as follows: 

-.290(CH1) - .562(CH2) + .600(CH3) + .491 (CH4) 

For any pixel, a greenness measure could be obtained, by using the above green- 
ness vector with the CHI - CH4 digital values of that pixel. This greenness 
measure is now called the Green Vegetation Index (GVI), and as seen from 
above, 

GVI » -.290(CH1) - .562(CH2) + .600(CH3) + .491 (CH4) 

Figure 2 clusters which were not soil but did not contain pixels having high 
greenness but instead contained what was called yellow stuff were selected, 
(actually fig. 2 did not contain any yellow pixels, and these had to be simul- 
ated from other scatter plots). Using Gram-Schmidt again, Kauth and Thomas 
created a new "yellow" vector which is orthogonal to both the soil brightness 
and greenness vectors. As a unit vector, this vector is described as follows: 

-.829(CH1) + .522(CH2) - ,039(CH3) - .194(CH4) 
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For purposes of finishing a matrix described later, and to describe a fourth 
factor Mhich cannot be defined as soil brightness, greenness, or yellowness, a 
fourth vector was created, and was said to be a feature of “none-such" or 
"none of the above factors". It Is singly a vector chosen to be orthogonal to 
the above 3 vectors, and was created by Gram-Schml dt again. As a unit vector, 
this "none-such" vector Is described as follows: 

.223(CH1) + .012(012) - .543(013) + .810(014) 


Using the four vectors above. It was now possible to define four Independent 
factor measures for any pixel, even though It Is accepted that only the first 
two vectors provide meaningful Information. Rather than using the vectors 
totally Independent of each other, Kauth's and Thomas' whole Ingenious system 
can be explained completely In terms of a single transformation, which Kauth 
and Thomas outline (ref. 13, page 48-43) as follows: 

Let u * R^X + r 

where x Is the Landsat MSS signal vector expressed In counts 

u Is the transformed vector, also expressed In counts 

r Is an arbitrary offset vector, simply introduced to avoid negative 
values In u 

R is a unitary matrix, whose columns are simply the aforementioned unit 
vectors of soil brightness, greenness, yellowness, and none-such 
respectively . 


In much of the literature, the whole transformation system Is explained 
multiplication of a rotation matrix by CHI - CH4_this Is, In essence. 


matrix, often cal 

led K, 

Is equal 

to 1 

i 

/ .433 

.633 

.586 

.264 


f- .290 

-.562 

.600 

,491 

IN • K * 1 

* 

1-.829 

.522 

-.039 

.194 


^ .223 

.013 

-.543 

.809 


as a 
the 
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and CHI - CH4 multiplied by the first roM equals the SOI, CHI - CH4 multiplied 
by the second row equals the 6VI, CHI - CH4 multiplied by the third row Is a 
yellow ntmber, and CHI - CH4 multiplied by the fourth row Is a none-»ich 
number. As mentioned before, the SBI and the GVI are the only measures of 
Interest In most cases. Figure 4 shows how SBI and GVI may be viewed 
graphically In three-dimensional space. 


Often In the literature, especially since LACIE, a different Kauth rotation 
matrix Is seen. The K rotation matrix shown above was the matrix derived by 
Kauth and Thomas using Landsat 1 data only. Landsat 2 and 3 have since been 
launched, and due to sensor calibratioi differences, modified K matrices nave 
had to be developed, although they serve the exact same purpose. When usiny 
Landsat 2 digital values, which currently Is the standard, the following 
rotation matrix Is used: 


.332 

.603 

.676 

.263\ 

-.283 

-.660 

.577 

.388 \ 

-.900 

.428 

.076 - 

.041 1 

r-016 

.131 

-.452 

.882/ 

values 

, the following 

rotation matrix is used 

.386 

.742 

.842 

.279\ 

- .329 

-.812 

.719 

.412 1 

-1.044 

.527 

.095 ■ 

-.043 / 

- .019 

.161 

-.563 

.93// 

y the 

first row equals 

the SBI, etc., just as 


before. Interestingly enough, the methods used In creating the rotation 
matrices for Landsat 2 and 3 data were not the exact same iwthods used by 
Kauth and Thomas In creating the rotation matrix for Landsat 1 data, although 
they could have been. Since the basic purpose of Part 11 of this paper Is to 
present the concepts and mathematics of the original developmental work done 
on various VINs, the differences In the methods that Kauth and Thomas used to 
create the Landsat 1 data rotation matrix and the methods used to create the 
rotation matrices used with Landsat 2 and 3 data will not be presented here, 
except very briefly below. The rows of the rotation matrix K for use with 
Landat 2 are transposed unit vectors characterizing Landsat 2 data rotation 
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Greenness 





Figure 4.- A graphical view of the orthogonal relationship 
between the SB! and the GVI (Source: ref. 10, p. 8). 
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Into Tassel led cap data space particularly oriented to suit Lambeck's XSTAK 
algorithm (ref. 12, page 717— this reference also provides other references 
Into the exact methodology used). Roms 1 through 4 of rotation matrix K for 
use with Landsat 2 are s1ti4)1y multiplied by 1.161, 1.230, 1.246, and 1.U62 
respectively (calibration difference adjustments) to obtain the four rows of 
the rotation matrix for use with Landsat 3 data (ref. 8. page 2). 

The most common use of the Tassel led Cap Transformation system Is a pixel 
scatter plot of GVI versus S8I, where GVI Is the vertical axis and SBI Is the 
horizontal axis, (ref. 7, page 3-5). Using this scatter plot, a horizontal 
line is drawn below which are approximately one percent of the pixels. This 
line or value Is called the soil line (getting rid of the bottom one percent 
Is thought to protect against swampy areas and low outliers). Subtracting the 
soil line value from the GVI yields a result called the Green Numoer, or 
sometin«s called the Kauth Vegetation Index (KVI). As a formula, the Green 
Number, or KVI, is represented as follows: 

KVI » GVI - (Soil Line Calculation) 

Any pixel having KVI = 0 should be bare soil, any pixel with KVI > G should 
have some green vegetation, and any pixel having KVI > 15 is considered to be 
highly green. 

In 1977, while studying drought, Thompson and Wehnianen (ref. 2U) created the 
Green Index Number (GIN). The GIN is not considered to be a VIN, since it is 
not a pixel measurement, but can be thought of as a "segment VIN". The 
creators defined the GIN to be the percent of pixels in a segment with KVI > 
15, (ref. 20, page 203), computed as follows: 

GIN « 


where 

N « number of pixels with KVI > 15, and 22,932 ■ number of pixels In a 
segment . 

The GIN may be said to be the percentage of a segment with a nearly full cover 
of green healthy vegetation. 
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2.5 OEVELOPHENT OF THE PV I /WO DVI 

Richardson and Wlagand (ref. 16) conducted research In 1977 trfth hopes of 
creating procedures that would account for soil background variations 
factor known to hamper Interpretation of ve^tatlon surface reflectance. A 
year earlier Kauth and Thomas determined that the data space distribution of 
soil reflectance variation In Landsat data Is confined to a line (In two- 
dimensional data space) or a plane (In three-dimensional data space), and that 
reflectance variation of developing vegetation grows perpendicularly out of 
the plane of soils. Using this Information, Richardson and Wiegand set about 
research that would lead to the creation of the perpendicular vegetation Index 
(PVI7, PVI6) and the difference vegetation Index (DVI). 

First, they obtained Landsat digital counts (In each of the four Landsat 
bands) from highly reflective bare soil, low reflectance bare soil, cloud 
tops, cloud shadows, and water on April 2, May 17, June 4, July 10, 

October 17, and December 10, 1975 In Hidalgo and Willacy counties, Texas (if 
interested, exact counts are on page 1543, ref. 16). Next, for each pairwise 
Landsat band combination, a linear regression was run for purposes of deter- 
mining Kauth's line of soil. A table displaying their results Is shown in 
table 2 (copied from ref. 16, page 1543). 

At this point It was determined that It was not necessary to continue their 
study using all of the MSS band pairwise con4>1 nations, and the following logic 
was used to eliminate all but band combinations (5, 6) and (5, 7); band com- 
binations (4, 5) and (6, 7) were eliminated because bands within the visible 
and Infrared are known to be highly Intercorrelated, band combinations (4, 6) 
and (4, 7) were eliminated because they had lower correlation coefficients and 
higher standard errors of estimate than did (5, 6) and (5, 7), and finally 
because band cond)1nat1ons (5, 6) and (5, 7) had been found useful In past 
studies. The linear equations for these two band combinations were designated 
as soil background lines, (Note: for (5, 7) the Intercept was not 

statistically significant; thus set to 0). 
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TABLE 2.- LINEAR EQUATIONS DETERMINING KAUTH'S LINE OF SOIL 
FOR ALL POSSIBLE PAIRWISE COMBINATIONS OF THE 4 LANOSAT 
MSS BANOS. DIGITAL COUNT DATA ARE FOR APRIL 2. mi 17. 


JUNE 4. JULY 10, OCTOBER 17. AND DECEMBER 10. 1975 
FROM HIGH AND LOW REFLECTANCE SOIL. AND CLOUD AND 
CLOUD SHADOWS (N - 16). 


MSS band 
pairwise 
combination 
(Xj. Xg) 

Correlation 

coefficient 

(r) 

Linear equations 
*1 " ®o * “ 1^2 

Standard error 
of estimate 
(Sxl-x 2 ) 

(4, 5) 

0.967 

Xj • -1.04 + 0 . 938 X 3 

Digital counts 
10 

(4, 6 ) 

0.949 

Xj - -5.45 + 1 . 011 X 3 

12 

(4, 7) 

0.958 

Xj • -1.23 + 2 . 257 X 3 

11 

(5. 6 ) 

0.993 

Xj • -5.49 f 1 . 091 X 3 

5 

(5. 7) 

0.987 

CM 

X 

0 

0 

CM 

0 

• 

0 

1 

II 

X 

6 

( 6 , 7) 

0.993 

Xj - 5.09 + 2 . 200 X 3 

4 


Richardson and Wiegand then plotted LAI values from some previously collected 
sorghum data for comparison with their soil bandground line (using the MSS5 
and MSS7 line). They noticed that data points for the sorghum fields deviated 
perpendicularly from the bare soil background line, and furthermore that the 
sorghum fields with larger LAI values (denser vegetation) wen displaced 
furthest from the line. They concluded that a measure of the distance of a 
candidate sorghtm point from the line could be used as an Index of vegetation 
amount for that sorghum point. They also noticed that water deviates fr'n the 
soil background line, but on the opposite side. Using the above findings, 
they concluded that Indeed the soil background line could perhaps serve as a 
soil background reference for a vegetation Index model. 

For both the (5. 7) and (5, 6 ) band combinations, a vegetation Index model, or 
VIN, was created simply using the perpendicular distance of a vegetation 
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candidite sl^atura point from tha soil oackground lino. Each will bo treated 
separately (to prevent confusion), using the notation and descriptions of tne 
originators. For MSSS Md HSS7, this perpendicular distance Is given by the 
equation: 

PVI7 -ll (R«„5 - R„5)^ ♦ (R„„7 - K„7)^ 

1 ' gg p ' ' gg p 

where 

PVI7_1s the perpendicular VIN, defined as the perpendicular distance 
between the candidate vegetation point and the soil background 
line for Landsat bands HSS5 and MSS7. 

Rp _ Is the reflectance of a candidate vegetation point for nsbb and 
M5S7, and 

Rgg _ Is the reflectance of soil background corresponding to a candidate 
vegetation point, 

and this distance shall be defined as “positive" If Rg^S > Kp5; “zero" if 
Rgg5 ■ Rp5, and "negative' If Rgg5 < Rp5. 

A picture displays this »«11: 



LANDSAT NSS7 
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since the toll tMckground line 1$ known. RpS end Rp7 for a candidate vegcta- 
t1c-«i point Is known (th^ are simply the Landsat digital counts for the vege- 
tation point), and the slope of the perpendicular line connecting (RpS. Rp7) 
with the soil background line Is known (It Is always the negative reciprocal 
of a line slope. In this case -1/2.4 ■ -.417), all that remains Is to solve 
for RggS and Rgg7, vdilch yields: 

R„5 ■ .851 R 5 ♦ .355 R 7, and 
99 P P 

■ .355 R 5 + .148 R 7. 

99 P P 

This, In turn, yields: 

PVI7 H ^'^99® ■ 

[(.850 Rp5 + .355 R^T) - R^S]^ + [(.355 Rp5 ♦ .148 Rp7) - Rp7)3^ 

- ^ (-.149 Rp5 + .355 Rp7)^ ♦ (.355 kp5 - .852 Rp7)^ 

■ ^ C.355(CH4) - .149(CH2)]^ + [.355(CH2) - .852(CH4)]^ 

which Is the form appearing in most of the literature. Using the above dis- 
tance formula, and the conclusions from the sorghum studies mentioned earlier, 
they were also able to conclude that a "positive" PVI7 Indicated vegetation, a 

"zero" PVI7 Indicated bare soil, and a "negative" PVI7 Indicated water. 

Richardson and Wiegand also created a con^)utat1ona11y slti^ler VIN using the 
soil background line of MSS5 and MSS7. It Is simply: 

DVI - 2.4MSS7 - MSS5 

- 2.4(CH4) - CH2 

and achieves the smne purpose as the PVI7. Just as with the PVI7, a "posi- 
tive" OVI Indicates vegetation, a "zero" DVI Indicates bare soil, and a "nega- 
tive" OVI Indicates water. However, with the OVI, the terms negative, zero, 
and positive have pure mathematical meaning: OVI < 0 Indicates water, etc. 
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How»>^r, ^ DVI hit a dlsadvantagi slfwa tht toll baclc 9 ^nd coordinatet 
(RggS, Rg^7) cannot bt determlnod. 

Anothtr ptrp«nd1cular VIN, PVI6, Mat alto crtated Alchardton and Wiegand. 

It Is computed exactly at PVI7 but uses Landsat bands 5 and 6 Instead of bands 
5 and K aid also uses the linear equation associated with bands 5 and 6, 

(HSSS • -5.49 4- 1.091 NSS6). This yielded: 


PVI8 -ij (H^i . Rp6)* ♦ (R„6 - Bp6)^ 


■ J C-.498 - .457(CH2) ♦ .498(CH3)]^ • [2.734 * ,498(CH2) - .543(CH3)]^ 

Richardson and Wiegand (ref. 16, page 1548) state reasons why they feel that 
the perpendicular VINs, PVI7 and PVI6, have a fundamental advantage over the 
other VINs existing at that time, which Included the TVI7, TVI6, OVI, GVI, 

SBI, and any single bands or simple ratio of two bands. Briefly, the PVI7 and 
PVI6 calculate the soil background Intersection coordinates which allow 
examination of reasons (water content differences, shadows, tillage, soil 
crusting) for differences In reflectance of cropland, rangeland, and forest 
scenes due to soil background. 

2.6 DEVELOPMENT Of THE LAI 

LAI, or leaf area Index, Is not a VIN but rather an Important agronomic 
variable. Evapotranspi ration ami photosynthesis models used regularly by 
researchers use LAI Inputs. Also, cr^p growth models use LAI Inputs 
Indirectly, since crop growth models typically Include evapotranspiratlon and 
photosynthesis subrmjtlnes (ref. 21, page 340). Correlation studies completed 
by AariHison and Davis (ref. 1) have Indicated a relationship between LAI and 
wheat yield. However, a large problem exists: grmjnd measurements of LAI are 
very tedious and expensive (and, of course, not possible In many foreign 
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artas). If LAI could cvtr bt estlMtad by usa of Landsat maasurtMonts , tna 
above mantlwtad usas of LAI could ba mada much mora Inaxpansivaly and could 
also ba axtandad to much largar araat. A fa« such afforts at LAI estimation 
have been made, almost exclusively by T. Kanemasu Morkiny «<1th nwstly «<heat 
dat4 


Standard regression analysis oas used to derive estimated LAI imidels In the 
knoMn efforts toward LAI estimation, vlth ground measurements of LAI as the 
dependent variable and Landsat band ratios or other VlNs as the Independent 
variables. Even though LAI Itself Is not a VIN, model estimated LAI nay be 
thought of as a VIN since It Is a number derived co«np1ete1y from Lanosat band 
combinations; thus model estimated LAI has appeared In swne VIN comparison 
studies . 

Although a number of LAI formulas have been proposed and used over the last 
few years, only three will be detailed In this section _ one because It was 
used In a large VIN comparison study (ref. 1) and at the time of this report 
Is still being used by FAS/CCAO In Houston, the others because they are 
currently accepted as the best wheat LAI estimates. All three are fo»‘ wheat 
only since this 1$ where the en^hasls has been. 

The LAI estimator used by FAS/CCAO during the last few years and by Aaronson 
and Oavis (ref. 1) in their VIN comparison study 1$: 

LAI . 41.325 (^) - 42.45 (^) . 

Unfortunately, this author has been unable to find a single reference with 
respect to the original development of this formula (Including its originator) 
or how well the model performed. It does appear similar to some of Kanemasu's 
earlier formulas; however, Kanemasu disclaims It. 
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In 1977, Kanemasu (ref. 21, page 339), using many observations taken over a 
tiKO year period from three large (over 40 ha) Kansas winter wheat fields, 
developed the following LAI model: 


( 1 , LAI • 2 . 67 ’ - 3.694 (^) - ..309 ® * 5.751 {^) 

+ nd? /CH2. , XH2 X . , n7i r^Hl CHI . .CHI. 

+ .043 (^) - 2.692 (^) + 3.071 (^ - ^) (^), 

An » .69 was obtained with this model. 


Kanemasu tried to improve the above model by using stepwise regression anal- 
ysis on 115 observations of LAI and their respective Landsat MSS data (again 
Kansas winter wheat) and a decision logic based upon low LAI versus high LAI 
areas (ref. 10, page 10). This resulted in the following model; 

CLAI = .366 - 2.265 (^) - .431 (^ - ^)(^) + 1.745 (^) .057 (PVI7); 

If CLAI < 0.5, then 


LAI * 1.093 - 1.138 (^) - .017 ( 


CHI 

CH2 


gi)(^) - .016 (PVI7), 


else 


LAI = -5.33 + .036 (PVI7) + 6.54 (TVI6). 

This model yielded an R^ = .69, no improvement over the of the simpler (1) 
model . 


Before leaving this section, it should be noted that Kanemasu has also 
produced various estimated LAI formulas for sorghum and soybean (ref. 11, 
page 46; ref. 21, page 338), however, these appear to be preliminary in nat’jre 
and not meant for usage. 
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3. VIN RELATIONSHIPS TO AGRONOMIC VARIABLES 


3.1 YIELD AND YIELD COMPONENTS 

This section discusses VIN relationships with yield and yield con^onents: 
percent crop cover, biomass, plant height, leaf area Index, and yield. 

Richardson and Wiegand (ref. 16) calculated correlation coefficients for the 
following Landsat bands and VINs with percent crop cover, plant height, and 
leaf area Index for ten sorghum fields: the Individual Landsat bands, band 5 
divided by band 7 (sometimes called RVI), TVI7, TVI6, PVI7, PVI6, DVI, GVI, 
and SBI. They found that the Individual Landsat bands had the highest corre- 
lations with all three agronomic variables. Band 5 was correlated highest 
with crop cover and plant height (r » -.809 and r » -.849 respectively), and 
band 6 was correlated highest with leaf area Index (r * -.849). The highest 
VIN correlations were TVI6 with crop cover (r * .716), TVI6 with plant height 
(r * .828) and PVI6 with leaf area Index (r = .812). All the above coeffi- 
cients were significant at the .01 level except TVI6 with crop cover, which 
was significant at the .05 level. 

Tappan (ref. 19, page 92) found that the single band ratios and band differ- 
ences seemed to correlate best with percent green vegetation cover when working 
with Kansas prairie data. All the single band ratios and s1n^)1e band differ- 
ences he worked with yielded r > .9. Tappan also noted (through personal 
communication) that percent ground cover and biomass have a direct relationship 
when vegetative ground cover Is less than 100%, and therefore the VINs which 
work well with percent cover also work well with biomass. Once 100% ground 
cover Is reached, band 7 has been found to be useful for detecting changes in 
biomass (ref. 19, page 25). 

Aaronson and Davis did a comparative VIN correlation study (ref. 1) using the 
AVI, DVI, GVI, KVI, LAI, PVI6, PVI7, TVI6, and TVI7, where LAI in this case Is 
the estimated LAI given by LAI * 41.325 (CH1/CH2) - 42.45 (CH1/CH3) described 
in section 8. Using winter wheat data, they correlated the VINs with each 
other, yield, plant height, and biomass at each of the following growth 
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stages: planting, tillering, stem extension, heading, flowering, ripening, 
ard harvest. Among their major conclusions were: 

1. All VINs In the study were highly correlated In each of the growth stages, 
(see ref. 1, page 10 for growth stage descriptions) 

2. The VINs were correlated to yield similarly at each of the growth stages, 
with highest correlations occurring at heading. LAI and TVI6 had the 
highest correlations with yield, and at heading had respective correlation 
coefficients of .64 and .63 (both at .0001 significance). 

3. The VINs were correlated to biomass and plant height similarly at tiller- 
ing and stem extension, with the highest correlations occurring at stem 
extension. Correlations with biomass at stem extension ranged from r = 

.65 to r ■ .71 (all at .0001 significance) for all the VINs with TVI6 
ranking highest. Correlations with plant height ranged from .76 to .84 
(at .0001) with AVI, DVI, and PVI7 ranking highest. 

Kanemasu, as noted in section 8, has been the pioneer in LAI estimation 
research, particularly for wheat. Some of his earlier estimated LAI models 
were singly functions of the PVI7 (ref. 21, page 337-338), and even his later 
two-step LAI model includes the PVI7 and TVI6 (ref. 10, page 10), thus showing 
the relationship noticed by Kanemasu between VINs and LAI. Working with 115 
observations of Kansas wheat data, Kanemasu found the following coefficients 
of determination between LAI and the below VINs and band ratios (ref. 10, 
page 11): 


VIN Band Ratio R^ 


PVI7 

.55 

CH2/CH3 

.56 

PVI6 

.55 

CH2/CH4 

.48 

TVI7 

.50 

CH1/CH4 

.27 

TVI6 

.59 

CH1/CH3 

.20 

GVI 

.57 

CH1/CH2 

.10 


3.2 STRESS FACTORS 

Three common stress factors affecting crop condition are water stress, 
nutrient deficiencies, and diseased or infested crops. The use of Landsat 
VINs In detecting these factors Is the topic of this section. 
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Thonpsofl and Uetmanen (ref. 20) successfully used the GIN in operation mode to 
detect and monitor drought In the U.S. Great Plains, U.S.S.R., and Australia 
during the 1977 crop year. In the U.S., they were able to statistically test 
their results against results using the Crop Moisture Index (CMI, ref. lA), 
which is known to detect water stress well, and found the drought test proce- 
dure based on the GIN worked very well. In the foreign areas, with no ground 
truth and CMI data available, they were unable to test their results statisti- 
cally, however, the U.S. agricultural attache for the U.S.S.R. verified that 
the U.S.S.R. regions which were Identified as water-stressed by the GIN test 
had Indeed undergone water stress at that time, and also a general moisture 
condition map produced by the Australian government showed ^neral agreement 
with GIN test results. 

Landsat D, when launched. Is scheduled to have thermal Infrared bands. It Is 
thought that these bands may prove to be useful in water stress studies since 
laboratory studies Indicate strong water absorption of Irradlance In the 1.2 - 
2.5 urn wavelengths (ref. 6, page 16-4; ref. 7, page 2-3). 

Research In the area of nutrient deficient crops has not yet reached the point 
of using any of the common Landsat VINs (not counting individual band VINs). 
However, many results have been obtained which are useful to know and cer- 
tainly will affect future Landsat research and desigi of new Landsat systems. 

Much of the study on crop nutrient deficiency has been done with corn. Al- 
Abbas et al. (ref. 2) conducted a study using nitrogen-, phosphorus-, potas- 
sium-, sulfur-, magnesium-, and caldum-deflclent corn leaves along with a 
control group of normal com leaves. He found that the chlorophyll concentra- 
tion In all nutrient-deficient leaves was lower than that of the normal 
leaves. He also found that in the near-infrared wavelengths leaves from the 
phophorus-and calcium deficient plants absorbed more energy than those frwn 
normal plants whereas leaves from the sulfur-, magnesium-, potassium-, and 
nitrogen-deficient plants absorbed less. 
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Bauer (ref. 5. pe^ 6} studied nltro^nt-deflclent com versus norMl com end 
found nitrogen deficiency caused Increased reflectance In the visible 
wavelengths and (tecreased reflectance In the near-infrared wavelengths. 

Cate, Artley, and Phlnney (ref. 7, pages 4-2, and 4-3) report that nitrogen 
deficiency has been f(Mind to Increase leaf reflectance In the visible bands in 
cotton, cabbage, and sweet peppers also. This same report also includes 
discussion on more nutrient (teficlency studies Involving mexican squash and 
sorghum. 

The use of Landsat In diseased or Infested crop research appears to be pretty 
dim for at least the Inmedlate future, a1thouc,h again there has been research 
which will affect future Landsat research and design. Studies conducted by 
Ausmus and Hllty (ref. 4) and Bauer (ref. 5, page 6} suggest that the near- 
infrared wavelengths (0.8 - 2.6 urn) were helpful in southern corn blight 
detection. Other studies (see ref. 7, page 4-5} also Indicate that near- 
infrared bands may be useful In disease detection. The launching of Landsat 0 
should help out in this area; however, Heller (ref. 9) advises that Landsat 
bands need to be both narrower and more selective in both the visible and 
near-infrared wavelengths to effectively detect vegetation damage, whether it 
be caused by disease or Insect-Infestation. A major "new" Landsat band sug- 
gested by Heller is a 0.58 - 0.62 band (a narrow yellow-orange band considered 
very useful for vegetation damage assessment). 

Disease and insect-infestation trends can be identified and monitored now by 
use of remote sensing, according to Heller (ref. 9, page 1159), bjt require a 
series of acquisitions over one or more seasons. 

3.3 CROP IDENTIFICATION AND GROWTH STAGE 

Using crop calendar knowledge, VIN computations considered over time have been 
used to help Identify general crop types such as spring small grains versus 
other, winter grains versus other, and summer crops versus other. Ashburn 
(ref. 3) discusses a nwjor well documented method which has been used in the 
past and also considers approaches to this problem in general. 
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Discrimination of specific crops within a general crop type (such as spring 
wheat versus spring barley) Is a further and sometimes difficult problem. For 
Instance, although spring wheat and spring barley have known differences In 
phenology when grown under Identical conditions. In actual cases these 
differences may te difficult to see due to differences In planting date, 
stress factors, and other factors. 

Any use of a VIN approach for general or specific crop identification depends, 
of course, on knowledge of growth stage. Aaronson and Davis (ref. 1, page 11), 
In their comparative study of nine different VINs, note that all VINs have a 
similar relationship to growth stage. The graph below (copied from ref. 1, 
page 12) depicts this general relationship. 
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4. SUMHARY 


This report was Intended to provide some Introduction, history, conceptual and 
mathematical principles, and uses of the more conmonly used Landsat vegetative 
Indices of today. Listed on the next page are the VINs considered In this 
report. 
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LANOSAT VINS 


AVI ■ 2(CH4) - CH2 
DVI - 2.4(CH4) - CH2 

GVI - -.283{CH1) - .660(CH2) + .577{CH3) + .388(CH4)* 

SBI - .332(CH1) + .603(CH2) + .676(CH3) + .263(CH4)* 

KVI - GVI - (Soil Line Calculation) 

LAI . 2.677 - 3.6M(^) - Z.309@ ♦ 5.751(^) . ,043@ 

9 CQ9/CH2 . ^ , n?wCHl CHI ..CHI. 

-2.692(^) + 3.071 - -^)(^) 

PVI6 C-.498 - .457(CH2) + .498(CH3)]^ + [2.734 + .498(CH2) - .543(CH3)]^ 

PVI7 - \ C.355(CH4) - .149(CH2)]^ + [.355{CH2) - .852(CH4)]^ 

Tvifi . ./ 'Cri3 - Ch2 ~ 

1 CH3 + CH2 

TVI7 - XCHy- CH2 — 

1 CH4 + CH2 

♦For standard Landsat 2 
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